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We study the effect of including flavor changing neutral currents (FCNC) in the analysis of the 
neutrino signal of a supernova burst. When we include the effect of the FCNC which are beyond the 
standard model (SM) in the study of the MSW resonant conversion, we obtain dramatic changes in 
the Am 2 — sin 2 29 probability contours for neutrino detection. 
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It is well known that there is considerable discrepancy between the observed solar neutrino fluxes and their theo- 
retical predictions based on the standard solar models. To solve the longstanding solar neutrino problem Q, the most 
natural scenario is to include the effect of matter on neutrino oscillations (the MSW effect) || . The SM interaction 
between the neutrinos and the electrons in solar matter would enhance the v e <-> z^ jT conversion and that could 
account for the apparent deficit in the solar neutrino flux. If indeed the MSW effect is the cause of the solar neutrino 
problem, a similar situation would exist for the neutrinos produced in supernova explosions, such as those detected 
from the 1987A supernova. The effect of the MSW effect in the light of the SM on the supernova neutrino signal has 
been studied extensively 

Recently, there has been a lot of discussions on possible new neutrino interactions beyond those of the SM. If these 
^vq i interactions do exist, they would have to be included in the study of neutrino propagation in matter. Interest in this 
' subject has heightened especially after the recent announcement of observing an excess of large-Q 2 e + p scattering 
events at HERA ||, with a possible signature for leptoquark formation. Also, possible Flavor Changing Neutral 
Currents (FCNC) and Flavor Diagonal Neutral Currents (FDNC) arise naturally in a number of extensions of the 
SM, such as the Left- Right Symmetric Models Q and Supersymmetric models with R-parity violation [Q. In fact, 
they have been recently discussed in the study of solar neutrinos p|-pd[| . 

f \ In this paper, we study the effect of these new physics interactions on the supernova neutrino signal. It will be 

Q\ • another test of these new models which covers a wider range of neutrino parameters compared to the solar neutrino 
""p^ | analysis. We neglect the non-standard FDNC interactions as the current limits on them make their effect negligible if 
neutrinos have mass. In section 2, we use the three-flavor formalism to derive the survival probability for neutrinos in 
the presence of FCNC beyond the SM. The three-flavor formalism is appropriate for describing supernova neutrinos 
Oh, because of the large range of densities in the supernova which the neutrinos must traverse. This means that, for all 
conceivable neutrino mass differences, both v e <-> and v e «-> v r conversions will go through the resonance regions ||. 
In section 3, we show the results of our calculations for the suppression (enhancement) ratio of the neutrino events 
due to supernova bursts which can be observed by the SuperKamiokandc fll2| ] and the Sudbury Neutrino Observatory 
(SNO) [|13| (already in the construction phase). The latter detector has the advantage of being sensitive to all neutrino 
flavors, not just VeS. We use typical values for the initial fluxes of the different neutrino flavors produced during a 
supernova collapse based on theoretical models describing supernova collapse and which are in a good agreement with 
the 1987A event. In our dm 2 — sin 2 9 contours, we emphasize the case where we have only v — d scattering, however, 
other new physics interactions can be easily included in the analysis. 
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II. NEUTRINO FCNC INTERACTIONS 



The only relevant SM interaction which is effective in the MSW effect is the standard VF-exchange charged current 
v e e — > v e e. However, in supersymmetric theories without R-parity and other non-SM models, new types of possible 
neutrino interactions with matter enter the picture. We can divide these interactions into two types: Flavor Changing 
Neutrino Currents (FCNC) and Flavor Diagonal Neutral Currents (FDNC). To be specific, we will present an analysis 
based on R-parity violating supersymmetric models. However, similar results as well as our parameterization of the 
effective interactions will follow had we chosen to use other theoretical models. In particular, we will be interested in 
the following lepton-number- violating contributions to the superpotential Q: 

KjkLiLjE^, (2.1a) 
KjkLiQM, (2.1b) 
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where L, Q, E c , D c are the usual lepton and quark SU(2) doublets and singlets, respectively, and i,j, k are generation 
indices. Because of the contraction of the SU(2) indices, the coefficients Xijk should be antisymmetric under the 
exchange of i and j. We will assume the coupling to be real to conserve CP. The superpotentials, Eqs. (2.1), lead to 
the Lagrangians 



Xi 



+ ele&i + (gfcmM - (t <- j)] + H.C. 



(2.2) 



and 



A' 



ijk 
e l L d R u L 



+ (4)*K) C 4 



H.C. 



(2.3) 



These interactions can give rise to v — d and v — e scattering, via squark and selectron exchange, respectively, if the 
couplings Xijk and Ay & are nonzero. The current limits on Xijk and A^ fe can be determined from charged-current 
universality constraints, the branching ratios of lepton-number violating interactions, etc. A detailed account of the 
limits on X^k and X'^ k is given in Jl4j . 

To obtain the neutrino survival probability, we start with the three-flavor propagation equation for neutrinos (see 
Ref. H), 





(2.4) 



where A = 2E«/2GfN b is the SM squared-induced mass, E is the neutrino energy, N e is the electron density in the 
supernova core, while U is the three-flavor mixing matrix given by: 



U = exp (iipXj) r exp (i(pX$) exp (iwAa)] 
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Xi are the Gell-Mann matrices while ip, ip and u> are the mixing angles corresponding to 



and v„ 



oscillations respectively. For simplicity, we will set the CP violating phase, 5, to zero in our calculation. The new 



physics parameters B and C describe the additional FCNC contributions due to v e 
respectively. They are given in terms of the couplings as 



and v e <-> v T interactions, 



B = 2EV2(G e N N e + G U N N U + G d N N d ) 



= 2EV2G F N e 



N 

2e u + e d + (e u + 2e d )j^ 



(2.6a) 



C = 2EV2(G'£N e + G'SN U + G'$N d ) 



2EV2G F N e 



i + K + e' d + «, + K) 



Nn 

N P . 



(2.6b) 



where Nf denotes the number density of the fermion type /, G* N and G^ are the effective couplings of the reaction 

v e f — > v^f and v e f — ► v r f, respectively; while = G n /Gf{G^/Gf)- Note that we assumed an equal number 

density of protons and electrons inside the supernova core. The order of the new physics corrections can be at most 
10~ 2 and this is due to the current bounds from rare decays, atomic parity violation, etc. The largest corrections 
would be due to transitions between the first and third generations p0| . The relevant parameters in terms of the b 
squark mass, mz, are: 
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They are currently bounded from above at the order of 10~ 2 . 

To find the effective mixing angles in matter, we diagonalize M . We can rotate the flavor basis by the factor 
exp(— iipAs) exp(— iip\i) and find the eigenvalues of the transformed mass matrix 

M 2 = exp(—iip\5) exp(— iipX^M 2 ex\){i^p\-[) exp(i(/3As) 

1 / S - AC 2u + 2 AC 2 + d AS 2ul + d 1 AS 2ip + d 2 \ 
= - AS 2UJ +d l E + AC 2uJ d 3 , (2.8) 

1 \ AS 2v + d 2 d 3 2(m 2 + AS 2 ) -d J 



where 



E = m\ + m 2 , A = m 2 — m 2 , 
do = —2S 2ip (BS^, + CC^), 
d\ = 2C v (BC^p — CS^), 
d 2 = 2C 2ip (BSif, + CC^), 
d a = 2S V (BC^ — CS^,). 

We are using the same notation of Ref. |3| with C v = cos ip, S v = sin ip, etc. To zeroth order in Asin2ip, the first 
matter-enhanced mixing angle due to the first resonance region is given by 

Asin2w + A . n . 

tan2 Wm = — 1 2.9 

—A cos z ip + A cos 2tt> 

with the new physics correction Di = d\ = 2(B cos cos i/' — C smifj cosp). The lower resonance occurs at A cos 2 ip — 
Acos2lj. 

We now turn to the higher resonance region, we rotate the flavor basis by exp(—itpXj), we obtain for M 2 : 




2(A + m 2 S 2 ) + AC% AS 2uJ C v + b (m 2 - K/2)S 2v + b 



AS 2u C v + b £ + AC 2uJ -AS v S 2uJ , (2.10) 

l + 2m 2 C 2 / 



(m 2 - A/2)S 2v + h -AS v S 2u AS 2 + 2m 2 C 2 



where A = E — AC 2ui bo — 2(BC^ — CS^) and b\ — 2(BS,p + CC^p). The second matter-enhanced mixing angle is 
(to zeroth order in A sin 2u/2ml) 

m? sin 2ip + Dh 

tan2 ^= J + 4cos2^ ^ 

where Dh = b\ = 2(Ccos?/> + Psin-0) and the higher resonance is at A = m 2 cos2(p. Note that our expressions 
for the mixing angles lead to the SM three-flavor solutions Jl5[ when Di^ — > 0. If we assume a mass hierarchy, 
where m 2 3> m 2 3> m 2 , we can separate the two resonance regions, each resembling a two-flavor resonance and the 
three-flavor survival probability can be written as 

P z (y e -» u e ) = P l (y e -» i^P^K -» i/ e ), (2.12) 

where p and Ph are given by fil| : 

Pi(v e — > ^e) = ^ + f ^ — 0(Acos 2 p — Acos2o;)P c (7,cj)^ cos 2w cos 2cj m (2.13a) 

and 

Pfc,(z/ e -> i> e ) = ^ + Q - ©(-4 - ?7i 2 ; cos2< ( 5)P c (7, tp) \ cos2ipcos2ip m . (2.136) 
The level-crossing probability, P c , is given by 

P c = | ( V2 (x+)\Mx-)) | = cxp [-| 7 ^W] ; (2-14) 
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7 being the adiabaticity parameter and F{9) is a function dependent on the density profile which reduces to unity for 
a linear density function. Note that the adiabaticity parameter, 7, is modified from its MSW form, it is given by Jll]] 

1 = 1sm |l + e o cot20| 2 , (2.15) 

where 7sm = 2Ecos20(dN fdx)/N — ' wnne S = (A = m| — m\, m|) and 9 = (ui, ip) for the lower and higher resonances 
respectively. We used the two-flavor result as both Pj and Ph have the two-flavor form with just a redefinition of 
variables. The parameter eo is given by 

N„ 

£0 = £e + 2e n + 6d + (e u + 2td) — — \resonance- (2.16) 

* * e 

III. SUPERNOVA NEUTRINOS 

We will now apply our results to neutrinos emitted in a supernova explosion. As neutrinos come out of the core 
of the supernova, they would traverse through a medium whose density decreases to zero from a maximum value 
of about 10 12 gm/cm 3 . This means that the resonance condition, A = A cos 29, will be satisfied for all conceivable 
neutrino masses. Thus, unlike the case of the solar neutrinos, supernova neutrinos will go through two resonance 
regions as they stream out from the core and, hence, the three neutrino formalism presented in the previous section 
is necessary to describe them. 

During a type II supernova burst, we have two types of neutrino emission processes, neutronization and thermal 
processes. The former is due to electron capture process, e" +p — > n + i/ e , while the latter occurs during the implosion 
of the supernova inner layers leading to the thermal production of neutrinos of all flavors in equal amounts. We will 
assume a Fermi-Dirac distribution for both the neutronization and thermal fluxes, F® n and F®*, with the following 
initial temperatures and luminosities | fl7| 

L* a = L%, Tl = Tj = 3MeV, 

L* = £* = Ll, Tl = Tl = 6MeV, 

L\ = 2L X , T™ = 3MeV, (3.1) 

where the superscripts t and n stand for thermal and neutronization, while the subscripts e, x and a denote the electron 
flavor, muon or tau flavors, and any flavor, respectively. Of course, there are uncertainties related to these values due 
to the rarity of supernova events and the complication of theoretical calculations and there are now more detailed 
numerical calculations on the time evolution of supernova. However, our approximate approach seems reasonable 
since there are very few neutrino events available at the present and the results shown are quite insensitive to the 
model dependent quantities for a supernova. In addition, we will need the density profile of the core and mantle after 
the core collapse, it is given approximately by Q 

p(r) = C/r 3 , lfT 5 < p < 10 12 g/cm 3 , (3.2) 

with C varying weakly with r over the range 1 < C/10 31 g < 15. We assume the ratio N e /N n ~ 0.4 G3]. 

Knowing the survival probability, P 3 (y e — > v e ), we can plot the ratio of the number of events in a detector with 
oscillations present to the number of events without oscillations. For a light-water Cerenkov detector there are two 
kinds of neutrino events: inverse beta decay which is isotropic and the v — e scattering which is extremely directional 
and occurs for all flavors. These two classes of events combined with the two supernova emission processes leads to 
three relevant ratios to calculate. We show our plots for directional events from neutronization, the relevant ratio is 
given by j| 

_ [ dE v F® n [P J eda{v e + e) + (1 - P) J eda{v x + ej\ 
SK ~ fdE v F^fed*("e + e) ' ^ 

Here J eda represents an integral over electron energies (from to £„) of the differential cross section times e, the 
energy-dependent detector efficiency; v x represents or v T . P is the survival probability of electron neutrinos given 
by the three-flavor expression, Eq. (2.12), which depends on four independent parameters U e 2, U e s, S and m^; so, 
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in order to make a 2-d plot, we will have to constrain two of the four parameters. In the case of a D2O detector, of 
which the SNO detector is an example, we will have the following important reactions with the deuterium nucleus in 
addition to the usual e — v scattering: 



and 



v{v) + d — >• n + p + v(u), (3.4a) 



v e (p e ) + d->p(n) +p(n) + e (e+). (3.46) 



Of course, the cross sections for these processes are very difficult to obtain analytically and we will use fit functions in 
our calculations to represent the cross sections, a(y + d), given in p0[ . The ratio of events with oscillations to those 
without will be equal to unity for the first type of event as it is equally sensitive to all flavors of neutrinos while the 
corresponding ratio for the second type is given by 

J dE„e [PFl + (1-P) Ff] o(y e + d) 
RsNO ~ fdE„eF°ta(v e + d) ' (3 ' 5) 

where again e is the efficiency of our detector. We assume in our calculations e to be the same for the SuperKamiokande 
and SNO detectors. The efficiency of detection depends on the electron energy as pl| 

e{E) = 1-exp [-((£-4.2)/4) 2 ] , (3.6) 

where E is measured in MeV. 

Figs. (2) and (3) show the contour plots for the SuperKamiokande detector and the SNO detector. Plots for 
different values of = G%/Gf : the ratio between the new physics coupling to the SM coupling, are shown. Note 
that we get a similar behavior in the case of the study of solar neutrinos (see Ref. p!l|). We use the following 
constraints: mi = 0, (7713/7712) = 250 and {\U e 2\ 2 /\U e 3\ 2 ) = 100. We neglect the earth regeneration effects in our plots. 
In Fig. (2), we have the expected SM behavior with the upper part of the contour corresponding to the adiabatic 
solution while the lower left side of the triangles correponding to the non-adiabatic solution. Note the two resonance 
regions correponding to the conversions between different flavors. The isocontours for the two experiments are the 
same since we assumed the same efficiency for the two detectors which is only an approximation. Fig. (3) shows the 
effect of including non-SM couplings. At = ±10~ 2 , the shape of the contours is significantly different from the SM 
case. We have a central region occuring at all values of the mixing angle where the event ratio is maximally different 
from 1.0 (0.2 in the case of SK and 2.3 in the case of SNO). Note that we have an enhancement in the number of 
observed events in the case of the SNO detector. At = 10~ 4 , we regain the adiabatic solution at large 8m 2 , also 
the central region divides into two regions. For negative couplings, we have an additional bump due to the divergence 
occuring at tan 29 = —2ed (leading to 7 = Jsm)- Note that for large values of the mixing angle, the effect of the new 
physics coupling becomes negligible and the large angle solution in the far right part of the plots becomes the same 
as the SM solution. Finally as we reach very small values of the couplings, ~ 10 -6 , we start to regain the SM 
solution. 

We should note that the parameters (Sm 2 , sin 2 29) are already bounded by other experiments: CHOOZ j22], Su- 
perKamiokande studies of solar |23| and atmospheric pij neutrinos, etc. The current bounds can be summarized 
as: 

• Small angle solution for solar neutrinos based on the MSW effect: 

3 x 10~ 6 eV 2 < Am 2 < 1 x 10~ 5 eV 2 , 

4 x 10~ 3 < sin 2 26> < 1.3 x 10~ 2 . (3.7a) 

• Large angle solution for solar neutrinos based on the MSW effect: 

6 x 10~ 6 eV 2 < Am 2 < 1 x 10~ 4 eV 2 , 

0.5 < sin 2 29 < 0.9. (3.76) 
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• Vacuum oscillation solution for solar neutrinos: 

5 x 10~ n eV 2 < Am 2 < 1 x l(T 10 eV 2 , 

0.7 < sin 2 26» < 1. (3.7c) 

• Solution for atmospheric neutrinos: 

3 x l(T 4 eV 2 < Am 2 < 7 x 10~ 3 eV 2 , 

0.8 < sin 2 26» < 1. (3.7d) 
IV. CONCLUSION 

Supernova explosions provide us with a very useful laboratory for testing the ideas of neutrino oscillations in matter 
as after their production in a supernova core, neutrinos travel through a wide range of densities not available anywhere 
else. In this paper, we include non-standard FCNC neutrino interactions into the treatment of neutrino propagation 
in a supernova. We may parametrize these effects as four-Fermi effective interactions (~ eGpvvff). Depending on 
the process, stringent bounds already exist on e, typically of the order of 10 -2 or less. However, it is found that the 
changes in the rates of neutrino signals are appreciable even for e as small as 10 -6 — 10~ 4 . We showed numerical 
results in connection with the SuperKamiokande and the SNO detectors. It should be noted that the dependence of e 
is intertwined with that of the neutrino mixing angles. Our results are therefore especially useful when combined with 
other measures which determine the bounds on mixing angles or when more definitive information on e are available 
from direct measurements [^|. In conclusion, observing neutrino signals from supernovae would be another test for 
non-standard models which include FCNC interactions, such as the R-parity violating SUSY models. It may impose 
new constraints on the non-standard neutrino couplings or, more interestingly, signal new experimental evidence for 
these models. This treatment is complementary to other studies of FCNC interactions in supernova where the bounds 
on these interactions are inferred from their effect on electron antineutrino spectrum and r-processes |p5| . 
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Fig. 1: The survival probability P 3 (h / e — » v e ) versus the neutrino energy E for the case of = (solid curve), 
e<j = 0.01 (dashed curve) and = —0.01 (dotted curve). 

Fig. 2: Contour plot for the SupcrKamiokandc and SNO detector for e u = e<j = e e = 0. 

Fig. 3: Contour plots for the Kamiokande and SNO detectors for = ±10~ 2 , ±10~ 4 and ±10~ 6 . 



